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We use Nambu-Goto numerical simulations to compute the cosmic microwave background (CMB) 
temperature anisotropies induced at arcminute angular scales by a network of cosmic strings in a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) expanding universe. We generate 84 statistically 
independent maps on a 7.2° field of view, which we use to derive basic statistical estimators such as 
the one-point distribution and two-point correlation functions. At high multipoles, the mean angular 
power spectrum of string-induced CMB temperature anisotropies can be described by a power law 
slowly decaying as £~ p , with p = 0.889 (+0.001, —0.090) (including only systematic errors). Such 
a behavior suggests that a nonvanishing string contribution to the overall CMB anisotropies may 
become the dominant source of fluctuations at small angular scales. We therefore discuss how well 
the temperature gradient magnitude operator can trace strings in the context of a typical arcminute 
diffraction-limited experiment. Including both the thermal and nonlinear kinetic Sunyaev-Zel'dovich 
effects, the Ostriker-Vishniac effect, and the currently favored adiabatic primary anisotropies, we 
find that, on such a map, strings should be "eye visible," with at least of order ten distinctive string 
features observable on a 7.2° gradient map, for tensions U down to GU ~ 2 x 10~ 7 (in Planck units). 
This suggests that, with upcoming experiments such as the Atacama Cosmology Telescope (ACT), 
optimal non-Gaussian, string-devoted statistical estimators applied to small-angle CMB temperature 
or gradient maps may put stringent constraints on a possible cosmic string contribution to the 
CMB anisotropies. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 



The idea that all elementary-particle forces could sim- 
ply be different manifestations of a single underlying force 
gained much strength after the successful elaboration of 
a unified renormalizable theory of electromagnetic and 
weak forces If the quest for the ultimate grand uni- 
fied group of symmetries G has so far been unsuccessful, 
the mere idea of unification leads to strong theoretical 
predictions. In particular, if the particle physics interac- 
tions were unified at high energy, the expansion of the 
Universe should have triggered spontaneous breakdown 
of their symmetries [2} . Kibble showed that the induced 
phase transitions may form stable topological defects [3| . 
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Cosmic strings are the linelike version of such primor- 
dial vacuum remnants. Their energy per unit length U 
is directly linked to the energy scale of the phase tran- 
sition during which they were formed. Moreover, once 
formed, cosmic strings are stable and should still be 
present nowadays. Although an inflationary era occur- 
ring after a string-forming phase transition would dilute 
the defects enough to render them unobservable, defects 
formation at the end of inflation is a generic feature of 
particle physics motivated scenarios [3]. Motivated by 
the coincidence that strings formed at the grand unified 
theory (GUT) energy scale would induce density fluctu- 
ations with amplitude close to the observed amplitude 
of galaxy fluctuations, they were once considered as se- 
rious candidates to explain structure formation in our 
Universe However, the COBE data, associated to 
large-scale structure observations, already suggested that 
adiabatic fluctuations seeded galaxy formation Q, and 
the most recent high precision measurements of the cos- 
mic microwave background (CMB) anisotropies leave lit- 
tle doubt that this conclusion is right 7] . Quite recently, 



2 



cosmic strings have also attracted the interest of the fun- 
damental string theory community: the embedding of 
inflation in string theory may indeed produce another 
class of stringlike cosmological objects, dubbed "cosmic 
superstrings," that could be fundamental or Dirichlct 
strings associated with extradimensional branes. Al- 
though of intrinsic Planckian-like energy density, their 
effective four-dimensional mass would be redshifted in 
presence of warped extra dimensions Q . 

These theoretical considerations motivate the study of 
string- induced gravitational effects, such as lensin g o r 
gravitational wave emission @, [H [Tl|, [H [H, [H, GJ,lll] • 
Using CMB and large-scale structure data, an upper limit 
to a potential cosmic defects contribution to the primor- 
dial inhomogeneity can be derived [I?], EH- The corre- 
sponding upper bound on the dimensionless energy scale 
GU (in Planck units) of local strings currently ranges 
from 2 x I CP 7 to 10~ 6 and is a weak function of the model 
used to describe the string network. With the notable ex- 
ception of Abelian string networks, whose induced CMB 
power spectrum has been recently obtained on the cur- 
rently observable angular scales many of the CMB 
analyses performed so far use analytical or semianalytical 
defect models which may or may not accurately mimic 
the cosmological evolution of a string network. These an- 
alytical approximations are often introduced to circum- 
vent the difficulties associated with the highly nonlinear 
evolution of cosmic strings in an expanding universe. The 
theoretical understanding of cosmic string evolution in 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse is still an active field of research which has led to 
the development of numerical simulations in corp orating 
all the defect dynamics [H HH, [H [H [H [H Hf [H HI . 
A drawback is that they have a limited dynamic range 
in redshift making it impossible to directly compute the 
evolution of a network of strings from their formation to 
the present time. 

These numerical limitations have motivated the search 
for simple signatures in the cosmological observables, 
such as straight temperature steps in the CMB tem- 
perature maps |29j. With the advent of new arcminutc 
CMB experiments, such as the Arcminute Microkclvin 
Imager (AMI) the South Pole Telescope (SPT) [U 
or the Atacama Cosmology Telescope (ACT) [32|, look- 
ing for strings directly into these maps may be promising. 
However, the efficiency of such searches would be greatly 
enhanced by the theoretical knowledge of the exact string 
patterns that may be imprinted in the microwave sky. 

The goal of this paper is to generate realistic sim- 
ulated CMB temperature maps that include the effect 
of cosmic strings on arcminute angular scales. For this 
purpose, we use high-resolution numerical simulations of 
Nambu-Goto strings based on an imp roved version of the 
Bennett and Bouchet code [22], [28|. As already men- 
tioned, numerical simulations cannot probe a very wide 
range of redshifts. But as far as small angular scales 
are concerned, we can use the approach introduced by 
Bouchet [33| and compute the strings' evolution starting 



at the last scattering surface. Our method allowed the 
generation of 84 statistically independent maps, which 
we used to extract some basic statistical properties as- 
sociated with the string patterns. Special attention has 
been paid to quantifying the systematic errors induced 
by the numerical initial conditions. We find that ap- 
plying the gradient magnitude operator to our simulated 
CMB temperature maps leads to the reconstruction of 
the string shapes on the past light cone, which suggests 
that it may be of some help for direct searches. We then 
add Gaussian perturbations of inflationary origin and the 
Sunyaev-Zel'dovich (SZ) effects to our simulated maps. 
Convolving the resulting maps with a diffraction-limited 
beam based on the specifications of the ACT experiment, 
we find that the gradient magnitude operator exhibits the 
strings' signature down to GU = 2 x 1CU 7 . While we de- 
fer the use of more sophisticated statistical estimators to 
a future paper, the clear non-Gaussian patterns showing 
up in the maps suggest that these estimators will signif- 
icantly improve the current constraints on GU. 

The paper is organized as follows. In Sec. [Ill we present 
our numerical simulations of cosmic string networks and 
the method used to simulate pure string-induced CMB 
temperature maps. Sec. [TIT] is devoted to a basic statis- 
tical analysis conducted on 84 simulated maps. We give 
the mean probability distribution function of the temper- 
ature fluctuations induced by strings and the correspond- 
ing mean angular power spectrum, with an estimation 
of the various associated errors. Some properties of the 
gradient magnitude of the string-induced temperature 
anisotropies are also discussed. In Sec. IIV1 we examine 
the observability of strings for an ACT-like experiment 
when the primary anisotropies, SZ effects, and the effect 
of beam smoothing are taken into account. The suitabil- 
ity of the gradient maps as strings tracer is reexamined 
in this context. We present our conclusions in Sec. El 



II. SIMULATED MAPS 

A variety of methods have been used to compute the 
CMB anisotropies induced by cosmic strings. Many ap- 
proaches rely on the use either of Green functions [3J] or 
of unequal time correlators (UTC) [HI], and all of them 
require the knowledge of the defect stress tensor evolu- 
tion during the cosmological expansion. Because of the 
intrinsic nonlinear evolution of defect networks, this is 
usually achieved through numerical simulations in FLRW 
space-time [20L |2l l l23l. [24l I2H . although simple analytical 
defect models [3g, |37|, |38|, [39( have widely been used to 
circumvent the numerical limitations mentioned in Sec. [J 

Computing the corresponding CMB angular power 
spectra requires only knowledge of the defect stress ten- 
sor two-point correlation functions. Since topological de- 
fects are active and incoherent sources of gravity pertur- 
bations, their two-point functions are generically nonva- 
nishing at unequal times |4(|. The UTC method relies 
on the scaling properties associated with an evolving cos- 
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mological string network to make this calculation easier. 
The network scaling properties indeed strongly restrict 
the UTC's functional form. They are then fit to numer- 
ical simulations and extrapolated to cosmological scales. 
This method has been successfully applied to derive the 
CMB and matter power spectra for a variety of topolog- 
ical defects [4l|, [42], EH, |3| and has recently been used 
with Abelian strings in Ref. 

Since we are interested in generating realistic maps, 
we need to go beyond computing the power spectrum. 
The two-point functions do not encode the non-Gaussian 
features induced by networks of defects in the CMB. We 
therefore use numerical simulations to capture the non- 
linear effects associated with the defect evolution, includ- 
ing non-Gaussianity. This approach has already been 
applied to constrain the energy breaking scale associated 
with various cosmic defects and to simulate full-sky CMB 
maps for Nambu-Goto strings [H, [H, S3, |H| . But due 
to the rather small expansion factor numerically achiev- 
able, such maps can only include stringy effects up to a 
finite redshift, typically z ~ 10 2 . The CMB anisotropics 
computed in this way are therefore only accurate on large 
angular scales. We avoid this limitation by stacking maps 
from different redshifts, an approach outlined in Ref. [33] 
and applied in Ref. [491]. While simulations with the ob- 
server outside of the numerical box are not well suited 
for a full-sky map reconstruction, they are quite useful 
for the small angular scales considered in this paper. 

Denoting Oi the £th multipole moment of the tempera- 
ture perturbation to photon distribution, the Boltzmann 
hierarchy in Fourier space can be recast into [50l . l5~fl ] 

no 2,2- 

H, = / fl(r?) e~ k / fc ° (e + <l) h(kAr]) dry 

lo 2 2 

r 1 g(ri)e~- k /k viv h j' e (kAr])dr] (1) 

n° _ T /d$ d*\ . A . 

where rj stands for the present conformal time, and $ 
are the Bardeen potentials, Arj = rj a — rj, and Vb refers 
to the baryon velocity A bar has been used when 

the Silk damping term has been explicitly extracted from 
the perturbation variables (53j. The visibility function g 
is related to the optical depth r through 



dr] 



(2) 



and is strongly peaked around ?7i ss , the time of last scat- 
tering. As a result, the contribution of the first two in- 
tegrals in Eq. |T]), namely, the Sachs- Wolfe and Doppler 
terms, is only significant for rj ~ 7yi ss [541 ]. Moreover, 
Silk damping exponentially washes these terms out for 
k > fcn, i.e., typically for £ > 2000. In presence of a cos- 
mological string network, the Bardeen potentials 'J and 
$ in Eq. ([T]) are sourced by all the usual cosmological flu- 
ids, as well as by the strings. Following Ref. 5], if \& s tg 



and <i> s tg refer to the perturbations that would be sourced 
by the strings alone, one can define = "J — ^ s tg 
and $ m = $ — $ s tg. As only a minor string contribu- 
tion appears to be compatible with the current CMB 
data [rEEl, ^stg < $stg < and one may consider 
the stringy effects as a first order correction to the stan- 
dard adiabatic cosmological perturbations. Although, in 
the presence of strings, the Bardeen potentials ^ m and 
<& m are coupled to \& and $ through the perturbed Ein- 
stein equations, at leading order, 
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and 
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(3) 



where the index "coh" refers to the purely coherent 
Bardeen potentials of inflationary origin obtained with- 
out strings. Therefore, dropping the Sachs-Wolfe and 
Doppler terms in Eq. (Q]), one gets 
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e e ~ / e- T (<f> coh + y coh )j e {kArj)drj 
fe»i J 



(4) 



e T (^ stg + ^ ste )j e (kAr ] )di 1 , 



where a dot denotes a derivative with respect to rj. The 
first term represents the integrated Sachs- Wolfe effect for 
the standard cosmological fluids and would vanish in a 
purely matter-dominated universe. However, due to the 
recent domination of the cosmological constant and the 
existence of radiation residuals at the surface of last scat- 
tering, this term contributes significantly to the large 
angular scales, but is unimportant at small ones [50j. 
Moreover, as the Universe is optically thick before last 
scattering, Eq. (0| becomes 



e 



stg 



-i k x-. 



dr;, 



(5) 



where x 7 (t7) = xiAr] is the photon path along the line 
of sight. At small angular scales, one may therefore 
expect the strings' signature in the CMB temperature 
fluctuations to be dominated by their integrated Sachs- 
Wolfe (ISW) effect from the last scattering surface. In 
the following, we use cosmic string numerical simulations 
and the so-called small-angle approximation to simulate 
CMB temperature maps according to Eq. ([5]). 



A. Small-angle approximation 

For Nambu-Goto strings, <& s t g and ^stg are solution of 
the perturbed Einstein equations sourced by the Nambu- 
Goto stress tensor. In the temporal gauge (X° = rj), with 
a = U / \/—g, it reads 



I 



da e X»X V - - X'^X'" } S 3 (x - X) , (6) 



where e 2 = X' 2 /(l — X 2 ), U is the string energy per 
unit length entering the definition of the Nambu-Goto 
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action, and a prime denotes a derivative with respect to 
the string world sheet spacelike coordinate a [55| . In the 
small-angle approximation, which is well suited for an- 
gles typically smaller than the Hubble angular size at the 
epoch of interest, Hindmarsh, Stebbins and Veeraragha- 
van showed that, in the case of Nambu-Goto strings, 
Eq. (O can be simplified to [5f| |57| 



e 



%-kS.GU 
l 2 



(u • X) e" 



-il-X 



' T eda. 



(7) 



Xfl: 



The wave vector 1 denotes the transverse component of k 
with respect to the line of sight n, whereas, in the tempo- 
ral gauge, u encodes the string stress tensor distortions 
of the photon temperature and reads 



u = X 



(n ■ X') ■ X / 
1 + n-X 



(8) 



As can be seen in Eq. ([7]) , only the strings that intercept 
the photon path can imprint their signature in the CMB 
temperature fluctuations. As a result, the knowledge of 
u, and therefore of the string trajectories X, is only re- 
quired on our past light cone. In the context of string 
numerical simulations, the trajectories of all the strings 
are computed for all times. Therefore, to compute u, 
one only needs to determine what parts of the string 
network intercept our past light cone at a given time. 
Equation ([7]) was also at the basis of the earliest string 
maps presented in Ref. [33| . However, our source terms 
in Eq. (El) are slightly different than the ones used in 
Refs. [33l . l58j | and include an additional longitudinal com- 
ponent with respect to the string trajectory (5(|. As dis- 
cussed in Ref. [57j , this component encodes a logarithmic 
correction to the temperature fluctuations induced by the 
string curvature and may be significant for wiggly strings. 



B. Cosmic string evolution 

Since Eq. §5§ involves an integration from the last scat- 
tering surface to the present time, we are left with simu- 
lating the cosmological evolution of a network of cosmic 
strings from z~1100toz~0. 

Our numerical simulations in FLRW space-time are 
based on an improved version of the Bennett and Bouchct 
Nambu-Goto cosmic string code [ID, [28j. The runs are 
performed in a comoving box with periodic boundary 
conditions and whose volume has been scaled to unity. 
The scale factor is initially normalized to unity, while 
the horizon size dh is a free parameter controlling the 
initial string energy within a horizon volume. During 
the computation, the comoving horizon size grows and 
the evolution is stopped before it fills the whole unit vol- 
ume to avoid spurious effects coming from the boundary 
conditions. We use the Vachaspati-Vilenkin initial con- 
ditions for which the long strings path is essentially a 
random walk of correlation length £ Cl together with a ran- 
dom transverse velocity component of root mean squared 



amplitude 0.1 [59j |. The initial horizon size and velocity 
amplitude are chosen to minimize the relaxation time of 
the Vachaspati-Vilenkin string network toward its stable 
cosmological configuration. As mentioned before, the fi- 
nite numerical box limits the range of scale factors during 
which the string evolution can be computed at once. In 
our case, to model the network evolution since the epoch 
of last scattering, we would need to follow it over a period 
during which the Universe expands by a factor of order 
10 3 , which is not achievable, even with the current su- 
percomputing facilities. Moreover, applying the results 
of the small-angle approximation would turn out to be 
difficult in such a case. Indeed, if z\ denotes the red- 
shift at which we start the numerical computation, the 
numerical box corresponds to a real comoving size 



2/Am. 



(9) 



in the matter-dominated era and for a flat universe. The 
density parameter Q m refers to the matter content of 
the Universe today, whereas dh is expressed in units of 
the numerical box size. For the largest simulations to 
date, z\ — 1089 and c?h — 5 x 10~ 2 are typical val- 
ues of these two parameters (see, e.g., Ref. |28|). In a 
flat ACDM universe, using fiducial values for the den- 
sity parameters compatible with the three-year Wilkin- 
son Microwave Anisotropy Probe (WMAP) data Q, this 
gives -L s im ~ 9 Gpc and the numerical box then subtends 
an angular size of 9f ov ~ 27°. The largest structures in 
such a map would therefore be out of the small-angle 
limit we are interested in. As a result, we adopt the ap- 
proach introduced in Ref. [H[ , which relies on two smaller 
runs. The first one starts at the last scattering surface 
and ends at a redshift fixed by the maximum expansion 
factor achievable in the numerical box. For the simula- 
tions we performed, dh — 0.185, which corresponds to 
L S i m ~ 1.7 Gpc and 8f ov ~ 7.2°. Such a run ends after 
a 30-fold increase in expansion factor, corresponding to 
a redshift z ~ 36. We then propagate the photons per- 
turbed by the first run into a second numerical simulation 
of the same size but starting at z\ ~ 36. For another 30- 
fold increase in expansion factor, this run ends at z ~ 0.3. 
As can be seen in Eq. ([9]), the second simulation repre- 
sents a much larger real volume than the first one and 
therefore subtends a greater angle in the sky. As a re- 
sult, only the subpart of the second run that matches the 
angle subtended by the first simulation is actually used. 
As we will see later on, the CMB temperature maps are 
weakly sensitive to the string network at low redshifts, 
simply because there are almost no strings intercepting 
our past light cone in a recent past, which makes this 
technique perfectly acceptable. 

In practice, each of these numerical simulations is 
started before the redshifts we mentioned in order to give 
the cosmic string network enough time to relax toward its 
stable cosmological configuration. As soon as this scal- 
ing regime is established, the string characteristic lengths 
evolve in a self-similar way with respect to the horizon 
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size, allowing the numerical results to be rescaled to cos- 
mological distances. Therefore, one has to make sure that 
the structures (strings and loops) we are interested in 
have indeed reached their scaling behavior during the nu- 
merical runs. This is not an issue for the so-called infinite 
(or long) strings, defined as strings larger than the hori- 
zon size, because they rapidly reach the scaling regime. 
Although it has been shown in Ref. [28| that the cosmic 
string loop distribution scales as well, the relaxation time 
for the loops to reach such a self-similar evolution with 
respect to the horizon size appears to be larger for smaller 
loops. As a result, and this is inherent to all cosmic string 
numerical simulations, the smaller length scales in a nu- 
merical string network keep some memory of the initial 
network configuration until they reach their stable cos- 
mological evolution (see also Re fs. [27|, |60|, l6l| ) . Note that 
even if this memory effect is physical, one does not expect 
a physical string network at the last scattering surface to 
still exhibit structures coming from its initial configura- 
tion at the GUT energy scale. The change in scale factor- 
between the GUT redshift and the last scattering surface 
is indeed so huge that all the observable length scales 
should be in scaling at decoupling. Unfortunately, as 
already mentioned, numerical simulations cannot probe 
the evolution of a cosmic string network up to the GUT 
epoch. To circumvent this issue, we switch on the photon 
propagation inside the runs only after making sure all the 
large structures (infinite strings and loops) are in their 
scaling regime. This can be checked, for instance, by 
monitoring the evolution of the energy density distribu- 
tions. In the following, we start the photons' propagation 
when all loops larger than a third of the horizon size are 
in scaling. This cutoff is then dynamically pushed toward 
smaller values to include all the loops entering the scaling 
regime at later times. The cutoff time dependence can 
be deduced from the loop distribution relaxation times 
derived in the simulations performed in Ref. [28]. We 
discuss the residual systematic errors associated with the 
presence of nonscaling structures in Sec. Ill CI Finally, let 
us point out that our Nambu-Goto simulations are not 
smoothing out the so-called "small-scale structure" along 
the strings [22j| . As a result, the CMB maps we generate 
depend on only one parameter, namely, the string energy 
per unit length U . 



C. Temperature fluctuation maps 

The CMB temperature fluctuations induced by strings 
on a 7.2° field of view are shown in Fig.[TJ The upper left 
(respectively, right) image corresponds to the tempera- 
ture fluctuations obtained at the end of the first (respec- 
tively, second) run mentioned above. The image in the 
top left panel exhibits more structures in the temperature 
patterns than the one in the top right corner. This prop- 
erty comes from the cosmological scaling of long strings 
which implies that their number inside a Hubble volume 
remains approximately constant during the expansion of 



the Universe. As a result, since the Hubble radius grows 
with cosmic time, the number of long strings intercepting 
our past light cone on a given field of view is higher at 
high redshift. It turns out that most of the visible long 
strings end up being located close to the last scatter- 
ing surface. Since this implies a small contribution from 
nearby strings, we did not perform a third simulation to 
fill the redshift gap between z = 0.3 and z — 0. 

The overall CMB temperature fluctuations are plot- 
ted in the bottom left panel, while the positions of all 
the strings intercepting our past light cone are repre- 
sented as a function of the redshift of interception in 
the bottom right corner. As can be seen on these im- 
ages, the temperature fluctuations appear as a superpo- 
sition of discontinuities associated with the long string 
segments. The presence of sharp edges is due to the 
motion of these segments with respect to the observer, 
which induces the typical Doppler shift associated with 
the Gott-Kaiser-Stebbins effect [9|, [lfj. Moreover, a few 
bright dipolar spots are visible and come from regions 
rich in loops and wiggly strings. These spots can be as- 
sociated with the cusps appearing on some isolated loops 
and with the kinks sourced in the active regions of string 
intersections and loop formation [5|[ . Although the over- 
all temperature map is dominated by the long strings 
pattern, the fluctuations induced by cusps and kinks may 
reach very high values. In fact, since the bright dipoles 
come from extremely localized regions, the amplitude of 
the temperature jump is smoothed on the scale of a pixel 
and increases with increasing image resolution until it 
reaches its physical value. As can be seen in Fig. [TJ the 
string-induced fluctuations observed with a resolution an- 
gle of 8 rcB = 0.42' (1024 pixels) clearly exhibit highly 
non-Gaussian structures. The basic statistical properties 
of these maps are further discussed in Sec. IIII1 



D. Robustness 

The CMB temperature fluctuations maps shown in 
Fig. [T] have been produced using the method described 
in Sec. Ill Bl to minimize the influence of the numerical 
initial conditions. The cutoff is adjusted at each redshift 
to keep all the loops which have entered their expected 
stable cosmological evolution and remove all the smaller 
loops which may come from the relaxation of the initial 
network configuration. A side effect of this procedure is 
to also remove the small loops which may be produced 
during the stable cosmological regime. Depending on the 
importance of these small loops on the resulting CMB 
temperature fluctuations, the simulated maps we obtain 
may therefore deviate from the ones that would be in- 
duced by a real string network. In order to quantify 
this deviation, we can compare our results (i.e., the map 
shown in the bottom left corner of Fig. [IJ thereafter, the 
"reference map" ) to the maps computed in two extreme 
cases: when all the loops are removed from the runs, 
and when all the loops, even the smallest ones believed 




FIG. 1: String-induced CMB temperature fluctuations on a 7.2° field with a (unrealistic) resolution of 8 lea = 0.42' (1024 pixels). 
The upper left image shows the fluctuations induced in between the last scattering surface and the redshift z — 36, while the 
upper right map represents the anisotropies produced by strings between z — 36 and z = 0.3. Because of their cosmological 
scaling, most of the long strings intercept our past light cone close to the last scattering surface. The overall string-induced 
fluctuations are plotted in the bottom left panel. As can be seen in the bottom right image, the edges in the temperature 
patterns of the other maps can be identified to strings intercepting our past light cone. Note that active regions corresponding 
to string intersection and loop formation events lead to the bright spots in these maps. Some of these spots are associated with 
O > 80 GU and saturate the color scale (see Sec. IIII|) . 



to come from the initial conditions, are kept. Those re- 
sults are shown in Fig. O As one may expect, the devia- 
tions are most apparent on the smallest scales. When all 
loops are included, the altered map differs from the ref- 
erence map by about 10%. The difference is a primarily 
small-scale noise with a few bright dipolar spots coming 
from the tiny loops. The situation is similar for the map 
obtained by removing all the loops. In that case, the 



few dipoles and large-scale structures produced by the 
loops in scaling are missing. Thanks to these "extreme" 
maps, we are now able to estimate the systematic errors 
associated with the presence of nonscaling structure in 
our numerical simulations. This is in particular useful to 
see how robust are the basic statistical properties of our 
maps, which we discuss in Sec. IIIII 
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FIG. 2: Influence of the nonscaling structures on the computed CMB temperature maps. a ii (respectively, Oi n f) refers to the 
CMB fluctuations that would be obtained by keeping (respectively, removing) all loops present in the numerical simulations, 
including the spurious ones coming from the relaxation of the numerical initial conditions. The left (respectively, right) panel 
shows the deviation induced by the presence (respectively, absence) of these structures with respect to the reference temperature 
map shown in the bottom left corner of Fig. [1] We use these maps in Sec. IHII to estimate the systematic errors induced by the 
presence of nonscaling structures in the cosmic string simulations. 



III. BASIC STATISTICAL PROPERTIES 

The CMB temperature patterns are directly related 
to the properties of the strings intercepting our past 
light cone in a particular realization of the associated 
network (see Fig. [1}. In order to extract meaningful 
properties from the string-induced CMB anisotropics, we 
generate 84 statistically independent maps following the 
techniques described in Sec. [Til Two CPU years have 
been devoted to performing 28 runs starting from statisti- 
cally independent Vachaspati-Vilenkin initial conditions 
(on IBM Power4 P660+ and AMD x86-64 2 GHz proces- 
sors). Each of these numerical simulations gives three 
temperature maps along the three spatial directions that 
can be used as either the low or high redshift contribu- 
tion. Finally, the overall temperature anisotropics map 
is obtained by combining a low and a high redshift map 
coming from two different runs. 



A. One-point distribution function 

The non-Gaussian nature of the temperature fluctua- 
tions induced by strings produces a very characteristic 
pattern on the maps shown in Fig. [1] A simple quanti- 
tative test of this feature is given by the one-point prob- 
ability distribution function (PDF) of the temperature 
anisotropics. The normalized PDF has been derived in 
three cases: temperature anisotropies generated by the 



infinite strings only (0i n f map), by the strings and loops 
in scaling, and by all the structures (O a ii map). The cor- 
responding PDF's averaged over all the simulated maps 
are shown in the left panel of Fig. [3] In each case, the 
PDF unambiguously deviates from a Gaussian as can 
be seen by comparison to the best Gaussian fit to our 
best PDF estimate (corresponding to the case where only 
strings and loops in scaling are taken into account). It is 
also worth pointing out that the result is relatively insen- 
sitive to the inclusion of loops, although including non- 
scaling structures increases the width of the distribution. 

More quantitatively, the mean of the distribution is 
zero to better than one part in 10 5 , whereas the variance 
averaged over all 84 maps is (ct s 2 ) ~ 154lg 4 (GU) 2 T| MB . 
The systematic errors quoted in this result have been 
derived from the mean values of the variance coming from 
the @ a n (upper bound) and 0i n f (lower bound) maps. 
A power-law fit of the tails of the PDF's restricted to 
fluctuations with |0| > 50 GU shows that they vary as 
0~ q with q = 8.2 ± 0.2, which, although steep, decays 
much more slowly than a Gaussian. 

As can be seen in the shape of the PDF compared 
to the best Gaussian fit, we have also measured a slight 
negative mean sample skewness 



(fJi 



(6 - 6) 3 



-0.24, 



(10) 



where (•) again refers to the average over all 84 maps, 
whereas a bar corresponds to an average over a single 
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FIG. 3: The left panel shows the probability distribution function of the CMB temperature fluctuations induced by cosmic 
strings (solid line). The dotted and dashed curves quantify the systematic errors coming from the string simulations. Each 
of these one-point functions is averaged over 84 independent realizations. The dash-dotted curve represents the best Gaussian 
fit. Deviations from Gaussianity are clearly apparent in the tails of the distribution, as well as from the slight skewness. The 
right panel shows the probability distribution function of the CMB temperature fluctuations that would have been induced by 
the nonscaling structures, as defined in Sec. Ill Dl Once again, the dash-dotted curve represents the best Gaussian fit. A slight 
positive skewness may be observed, suggesting that the negative skewness observed in the left panel is the result of the loop 
formation mechanism. 



map. Although the dispersion of g± between the maps 
is significant, there is a clear preference for g± < with 
less than 15% of our maps exhibiting a positive skewness. 
This skewness may come from high correlations between 
the string velocities and the shape of the strings along 
the past light cone. To confirm this, we checked that the 
skewness disappears if the velocity field along the long 
strings is randomly redistributed. Since the shape and 
velocity field of strings strongly depend on their non- 
linear evolution, which itself relies on intercommutation 
and loops formation, one may interpret this effect as a 
consequence of their nontrivial interactions. This inter- 
pretation can be further explored by looking at the PDF 
of the temperature fluctuations associated with the "non- 
scaling" structures, as defined in Sec. Ill Dl This PDF is 
shown in the right panel of Fig. [3] and exhibits a slightly 
positive skewness. Two kinds of structures contribute 
to this PDF: the tiny loops generated during the sta- 
ble cosmological evolution of the network, and the loops 
generated by relaxation of the initial conditions. It is 
precisely because of the difficulty of distinguishing the 
physical loops from the ones associated with the relax- 
ation of the initial conditions that we discarded these 
structures when generating our maps. As a result, we 
stress that using this PDF as a tool to understand the 
physics is not straightforward. The fact that the PDF 
shown in the right panel of Fig. [3] is also skewed sup- 
ports the existence of the correlations mentioned above. 
However, since the nonscaling structures are by definition 
still sensitive to the numerical initial conditions, the bare 



value of the skewness given in Eq. (fTOjl should be taken 
with a grain of salt. In particular, one may not exclude a 
smaller value if all the loops were evolving in their scal- 
ing regime, as it should be for a physical string network. 
To complement this qualitative explanation, it would be 
interesting to try to understand this effect quantitatively 
using analytical methods, such as the velocity dependent 
one-scale model of Refs. [62|, |63J. Let us finish this dis- 
cussion by pointing out that the value of the kurtosis for 
the PDF presented in the left panel of Fig. [3] is 

(g 2 )^ (e ~ 4 e)4 ^-3^0.63. (11) 

Gaussian fluctuations, such as the primary anisotropics 
of inflationary origin, are entirely described by their one- 
point distribution and two-point correlation functions. 
As the string-induced anisotropies are not Gaussian, the 
angular power spectrum does not contain all the statis- 
tical information associated with them, but provides a 
way of estimating their relative contribution to the CMB 
anisotropies as a function of angular scale. 



B. Power spectrum 

The two-dimensional power spectrum of the string- 
induced temperature anisotropies can easily be derived 
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FIG. 4: Mean (over our 84 maps) of the angular power spectrum of string-induced CMB anisotropics (solid lines). In the left 
panel, the dashed (respectively, dotted) line is the mean power spectrum obtained from string simulations including all loops 
(respectively, infinite strings only). These two curves give an estimate of the systematic error in the results derived from our 
numerical simulations. In the right panel, the dotted lines define the ±1-<t statistical error envelope around our best estimate 
of the power spectrum. The shaded area includes the effect of our systematic errors in addition to the statistical ones. 



in Fourier space. Assuming statistical isotropy, 

2 



P(u) 



I Bin. '«",!-« 



(12) 



where 6 = (a, (3) denotes angular coordinates and u the 
two-dimensional wave vector. In the small-angle approx- 
imation, this is also the angular pow er spectrum up to a 
redefinition of the wave number I64| . We have 



1(1+1) 
4tt 2 



C e ~ u 2 P{u), 



(13) 



with u = £/(2tt), which we use to produce Fig. 3J For a 
field of view of 7.2° and a resolution of 0.42', the multi- 
pole moments run from £ ~ 10 2 to 10 4 . As detailed in 
Sec. HIl we only take the string-induced ISW effect into 
account when computing the CMB temperature fluctu- 
ations from strings. As a result, the spectra we present 
in this section are only accurate at small angular scales, 
typically for multipoles higher than a few hundred. The 
left panel of Fig.[4]shows three spectra averaged over our 
84 maps. The central one corresponds to our best es- 
timate, whereas the top and bottom spectra have been, 
respectively, obtained from the string-induced tempera- 
ture anisotropies a n and 0j n f and trace the systematic 
error in the power spectrum estimation (see Sec. HID]) . In 
the right panel, we show our best estimate of the power 
spectrum, as well as of the associated systematic and sta- 
tistical errors combined. We define this combination as 
the smallest area around our best estimate containing all 
the ±1-(T statistical error envelopes of the spectra pre- 
sented in the left panel. 



Although the uncertainties induced by the presence of 
nonscaling structures in the simulations affect the overall 
normalization of the spectrum, its global shape remains 
the same. Using a power law fit to the small-angle tail of 
the spectrum, 



£(£+l)C e oc r p with 



p= 0. 



+0.001 
-0.090 ' 



(14) 



where we only quote systematic errors estimated from 
the extreme power spectra shown in the left panel of 
Fig. |4j Let us point out the fact that the overall am- 
plitude, in units of GU, of our best estimate spectrum 
appears to be compatible with the one recently ob- 
tained for Abelian strings at large angular scales by the 
UTC method 19]. At multipole £ = 1000, we have 
£{£ + l)C £ /(27r) ~ 14 {GU) 2 . It is nice to see these 
two techniques agree over the range of multipoles where 
they are both valid approximations, but it is not triv- 
ial that this should be the case. In particular, numeri- 
cal simulations of Abelian string networks do not exhibit 
the typical loop formation events observed for Nambu- 
Goto strings [3, H^]. However, the scaling regime of 
long strings is reached in both Abelian and Nambu-Goto 
simulations. As already discussed, the typical distance 
between long strings in scaling at the surface of last scat- 
tering is a fraction, typically a fifth [28), of the Hubble 
radius at that time. As a result, even the long strings 
are contributing to the CMB anisotropies down to mod- 
erately small angles. It is therefore not totally surpris- 
ing that both spectra are similar at intermediate angular 
scales (1000 < £ < 3000). We do not try to compare 
our results to those presented in [li| at larger angular 
scales since our method relies on the flat-sky approxima- 
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tion, which breaks down for multipoles smaller than a 
few hundred. 

As can be seen in Fig. [4J some extra power shows up 
at very small scales for the Cg's coming from the numeri- 
cal simulations including the nonscaling loops (0 a ii). The 
same effect was already visible on the corresponding tem- 
perature map plotted in Fig. [2] (left panel) and suggests 
that nonscaling structures start to have significant effects 
at very small scales, for t > 10 4 . The very small-scale 
effect of loops (be they in scaling or not) is also what ex- 
plains the asymmetric errors appearing in Eq. (|14p . As 
expected from the presence of cosmic strings all along 
our past light cone, the power spectrum decays quite 
slowly at small angular scales, and much slower than 
the exponential Silk damping experienced by the CMB 
anisotropics of inflationary origin. As a result, a contri- 
bution of cosmic strings to the angular power spectrum 
could be negligible at low multipoles and yet dominate 
the primary anisotropies for large values of I. We discuss 
the observability of string-induced anisotropies at small 
scales in more detail in Sec. IIVI 



C. Gradient magnitude as a string tracer 

Although the one- and two-point functions already ex- 
hibit some characteristic features of string-induced CMB 
anisotropies, they remain weakly sensitive to most of 
the highly non-Gaussian features one can see in Fig. [T] 
More sophisticated statistical tools have been developed 
to look for patterns in the CMB. These techniques in- 
volve the use of Minkowski functionals, the skeleton, 
wavelet analysis, or temperature step finders, among oth- 
ers [1^, [6|1 HI, [6?], The corresponding estimators 
are usually applied directly on CMB temperature maps. 
However, the deviations from Gaussianity appearing in 
the PDF shown in the left panel of Fig. are due to two 
main effects, namely, the presence of rare high tempera- 
ture events coming from localized string and loop regions 
with kinks and cusps, and the linelike (and fractal) shape 
of the temperature steps induced by moving strings in the 
microwave sky, which also explains the power-law behav- 
ior of the power spectrum at small angular scales. One 
may therefore expect such estimators to be more efficient 
when directly applied to a string map, such as the one 
shown in the bottom right panel of Fig.[TJ Of course, such 
a map cannot be directly observed. But since strings 
appear as discontinuities in the temperature maps, the 
temperature gradient would be singular at the string lo- 
cation, enabling us to produce a map similar to the one 
shown in the bottom right panel of Fig. [TJ though with no 
redshift information. This map could then be used with 
the techniques mentioned above to get more statistical 
information on the properties of networks of strings. 

Directional gradients were first discussed in Ref. [69[ 
to analyze the topology of string-induced temperature 
anisotropies. In order to conserve isotropy, it is con- 
venient to consider the gradient magnitude |V0| of the 
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FIG. 5: Normalized gradient magnitude of the string-induced 
temperature anisotropies shown in Fig. [T] (bottom left panel). 
A logarithmic scale has been used to enhance the contrast by 
preventing the bright spots from saturating the color scale. 
Such a map reproduces the string path on our past light cone 
and enhances the active string regions (see Fig. [TJ. 



temperature anisotropies defined by 

™-i/(£)'+(i)'. <-> 

where a and f3 are the horizontal and vertical angular co- 
ordinates. This definition makes it clear that for a finite 
temperature step, let us say 0(cv, (3) — © H(a; — a ), 
H being the Heaviside function, the resulting gradient 
magnitude is a Dirac distribution at the string location. 
On a pixelized map, the maximal amplitude of |V0| is 
therefore given by the size of each pixel. To take this 
resolution effect into account, we instead look at the nor- 
malized gradient magnitude |V n 0|/Gf , defined by 

|V n 0| = |V0|0 rcs , (16) 

of the temperature map shown in the bottom left panel 
of Fig. Q] The resulting map is shown in Fig. \5[ Note 
that we enhanced the contrast by using a logarithmic 
scale to prevent the bright dipoles from saturating the 
color scale. As can be checked by comparing Fig. [5] to 
the bottom right panel of Fig.[T] the gradient map repro- 
duces the string path on our past light cone. In addition, 
the magnitude of the gradient, once normalized to the 
image resolution, encodes the transverse string velocity, 
which renders kinks and cusps clearly visible. But since 
all the results presented in this section have been derived 
for strings alone, and with a high and quite unrealistic 
angular resolution, one may question their relevance in 
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view of the current experiments and data. In the next 
section, we address this issue by combining the string- 
induced anisotropies with the standard primary and ex- 
pected secondary CMB signals for a typical arcminutc- 
resolution CMB experiment. 



IV. OBSERVING COSMIC STRINGS 

CMB experiments have so far been able to map the mi- 
crowave sky down to scales corresponding to multipolcs 
of order 2 x 10 3 , WMAP providing cosmic variance lim- 
ited measurements up to I ~ 400 Q. At these scales, 
the primary CMB dominates over all secondary effects 
by more than an order of magnitude. However, upcom- 
ing experiments will produce CMB maps with arcminute 
resolution, at which point secondary effects become the 
dominant source of fluctuations [13, HH, H2] • It is there- 
fore crucial to take them into account to produce reason- 
ably realistic maps of what the CMB would look like at 
small angular scales in the presence of strings. Through- 
out this section, we consider a flat power-law ACDM cos- 
mology with h = 0.70, Q A = 0.73, fl h = 0.045, n 3 = 0.95, 
tt m = 0.27, cr 8 = 0.78 and r = 0.073. The hydrogen frac- 
tion is set to 0.76. 

A. Contribution of secondary anisotropies 

In the flat-sky limit, provided that the angular scale be 
small enough for the two-halo term to be negligible, the 
thermal Sunyaev-Zel'dovich effect (tSZ) angular power 
spectrum can be computed in the context of the halo 
formalism and is then given by [7(J [7l|, [73] 



C> 



1 



dz 



dV(z) 
dz 



M„ 



(17) 



where x — hv /k^T CMS , V{z) is the comoving volume at 
redshift z per steradian, dn/dM is the comoving number 
density of dark matter halos of mass M at redshift z, and 
yi(M, z) is the Fourier transform of the projected Comp- 
ton y parameter. To compute dn(M, z)/dM, we use the 
Jenkins mass function [731 ] . and the linear matter power 
spectrum produced by CAMB [74| . y&{M, z) is calculated 
as in Ref. [72|, the tSZ profile being integrated over a 
maximum of two virial radii. At the angular scales of 
interest, the values of z max , M mm , and M max have been 
set according to Ref. [72|] and are accurate enough to en- 
sure the convergence of the integrals. However, there are 
several theoretical uncertainties in the calculation of the 
mass function and the Compton y parameter, leading to 
a typical uncertainty of a factor of 2 on the overall am- 
plitude of the tSZ spectrum. Moreover, as the two-halo 
term is neglected, the current approach holds only for 
£ > Lax — 300 75]. This limitation is unimportant for 



our study as the primary CMB dominates over all other 
secondary effects by more than 3 orders of magnitude for 
multipoles smaller than £ max . 

The kinetic Sunyaev-Zel'dovich effect (kSZ) has been 
computed exactly in the linear regime, where it is also 
known as the Ostriker-Vishniac (OV) effect [zSlzSlzi. 
Although no exact result exists in the nonlinear regime, 
several analytical models have been proposed to describe 
this component [79|, H(| [1H- Here, we use the model of 
Ref. [8lj with, as input, the same linear matter power 
spectrum generated by CAMB as for the tSZ. Although 
the nonlinear corrections are unimportant at large scales, 
they become significant in the small-angle limit. Within 
our cosmological model, assuming instantaneous reion- 
ization, the nonlinear component of the kSZ has the same 
amplitude as the linear one at £ ~ 2500. 

Several other secondary effects can produce fluctua- 
tions in the CMB temperature, among which are gravi- 
tational lensing, the ISW and Rees-Sciama effects linked 
to the standard cosmological fluids, patchy reioniza- 
tion, and point sources. The corrections to the pri- 
mary anisotropies due to gravitational lensing have been 
included and computed with CAMB. Though they be- 
come important relative to the primary anisotropies for 
I > 2000, the Sunyaev-Zel'dovich effects are always dom- 
inant at these scales. The ISW and Rees-Sciama effects 
linked to the standard cosmological fluids are never dom- 
inant secondary effects for I > 100 [82j . and we therefore 
ignore them. It has been shown that patchy reionization 
can make a significant contribution to the nonlinear kSZ, 
though an agreement on how big this contribution is is 
still to be found [H, HH HH • Given the uncertainties in- 
volved in this effect and the fact that, even if the signal 
coming from patchy reionization dominates the nonlin- 
ear kSZ, it is still at least an order of magnitude smaller 
than the tSZ at all scales [85{ . we do not include it in our 
analysis. Finally, the question of how one should deal 
with point sources is probably the trickiest one. Mod- 
eling point sources (be they radio or infrared ones) is 
indeed an arduous task leading to quite uncertain mod- 
els [86| , which leads us to omit them in the present anal- 
ysis. This is a clear limitation of our study that should 
be kept in mind when trying to use the results presented 
here to predict how well we could be able to do as far 
as detecting cosmic strings is concerned. Although com- 
plementary observations at other wavelengths will help in 
removing the signal from the brightest point sources [321 ] , 
our ability to remove the remaining contamination will 
limit the possibility of detecting the signature of cosmic 
strings in the microwave sky. This should be yet another 
motivation for trying to improve our understanding of 
point sources and develop efficient techniques to remove 
them without losing interesting string-induced features. 

The angular power spectra of all these effects, namely, 
lensed primary anisotropies, tSZ, OV, and nonlinear kSZ, 
are shown in the left panel of Fig. [SI along with vari- 
ous spectra of the anisotropies induced by cosmic strings. 
From top to bottom, the string power spectra have been 
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FIG. 6: (Left) Angular power spectrum of the primary CMB anisotropies (solid line), the lensed CMB (dotted line), the thermal 
Sunyaev-Zel'dovich effect (long dashed-short dashed line), the Ostriker-Vishniac effect (lower long dash-dotted line), and the 
nonlinear kinetic Sunyaev-Zel'dovich effect (upper long dash-dotted line) in a flat power-law ACDM universe (see text). The 
fluctuations generated by cosmic strings are shown for three values of the string tension: GU = 2 x 10 -6 (short dashed line), 
7 x 10 -7 (long dashed line), and 7 x 10 -8 (short dash-dotted line). (Right) Model of the primary ACT beam in multipole space 
at 279 GHz (solid line), 215 GHz (dashed line) and 147 GHz (dotted line). The normalization is chosen such that the beam is 
unity at I = 0. As expected for a diffraction-limited beam, the cutoff multipole increases with increasing frequency. 



generated with GU = 2 x 1(T 6 , 7 x 1(T 7 , and 7 x 1(T 8 . 
The tension 2 x 10~ 6 is the one required to normalize the 
large-scale string-induced power spectrum presented in 
Ref. [H[ to the WMAP three-year data at I = 10. Such 
a high value of the string tension is already ruled out at 
more than 3-er [87] , but is commonly used as a reference. 
GU = 7 x lCT 7 is the 95% confidence level upper limit 
set in Ref. [8c| with the power spectrum presented in 
Ref. [l|| . Since the amplitude of our string spectrum at 
I ~ 1000 matches the one derived in this reference (see 
Sec. IIII Bp . we use GU = 7 x 10~ 7 to produce the maps 
shown in Figs. [5] and O Let us point out the fact that the 
associated string spectrum becomes marginally dominant 
over all primary and secondary anisotropies for multi- 
poles larger than typically 3 x 10 3 . However, for smaller 
tensions, the dependence of the amplitude of the string- 
induced power spectrum on (GU) 2 rapidly brings the lat- 
ter below the spectra of all the secondary anisotropies we 
consider for I above a few hundred. This is the case, e.g., 
for GU — 7 x 10~ 8 . Because of this strong dependence 
on GU, the factor of 2 uncertainty on the amplitude of 
the tSZ is relatively unimportant for our study. 

Combining all these results, we can now generate sim- 
ulated maps of the temperature fluctuations produced 
by the corresponding primary and secondary effects and 
add them to the string-induced anisotropies. However, 
to more realistically predict what arcminute-resolution 
CMB observations would look like in the presence of 
strings, we also need to convolve the resulting maps with 
a model for the beam of such an experiment. 



B. Instrument beam 

Several arcminute-resolution CMB experiments are al- 
ready underway, covering frequencies from 15 GHz (3fjj 
all the way to 345 GHz [3l| . In the remaining of this pa- 
per, we will use a model based on the ACT specifications 
to include realistic resolution effects in our maps. ACT 
is a 6 m telescope that will map the microwave sky at 
147, 215, and 279 GHz [II]. The actual data collected 
in such an experiment is not directly the temperature on 
the sky, but the result V of the convolution of the latter 
with the instrument beam such that 16411 



V(k) 



dB v {T) 
dT 



9(r)A(r) e" ik r d 2 r, 



(18) 



where B V (T) is the Planck function, r the coordinates 
of a point on the telescope, and A(r) the primary beam 
that we model as an Airy pattern. We are interested in 
the Fourier transform of the primary beam per unit area, 



A(u) = A 



u u u 

arccos \ / 1 — — 

u c u c V \ u c 



(19) 



where u = k/(2ir), A = 2/(n 4 d 2 ), and d = 6m. The 
beam in Eq. (j 19[) is defined for u < u c , where, at a given 
wavelength A, u c = 0/X, with 9 the characteristic maxi- 
mum opening of the telescope set to 70° (89[. For u > u c , 
A(u) — 0. The normalization is chosen such that, with 
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FIG. 7: CMB temperature anisotropics, including the lensed primary fluctuations predicted by our fiducial ACDM model (see 
text), as well as the tSZ, OV, and nonlinear kSZ effects, for a 7.2° field of view at 147 GHz (left), and corresponding normalized 
temperature gradient magnitude map (right). These maps take into account the effect of the primary ACT beam as derived in 
Sec. IIVBI The effect of strings is not included. 



the convention 



A{r) = jhf J A (u)e 2i<u r) d 2 u, 



(20) 



A(0) = 1. In the small-angle approximation, Eq. JT 
becomes, in multipole space, 



A(£) = A 



(21) 



with £ c = 2Ttd/{Xff). The normalized primary beam 
A(£) = A(£)/A(0) is shown in the right panel of Fig. © 
for the three ACT frequency bands. 

C. Generating mixed maps 

In this section, we present small-angle maps of the 
CMB temperature anisotropics induced by all of the ef- 
fects mentioned above. Producing these maps turns out 
to be relatively straightforward for two reasons: the pri- 
mary anisotropics of inflationary origin, but also the sec- 
ondary ones, namely, the tSZ and kSZ contributions, can 
be treated as Gaussian effects, and these fluctuations are 
uncorrelated with the string-induced anisotropics. The 
first of these two properties allows us to easily generate 
maps of the corresponding effects. Rocha et al. indeed 
showed in Ref. [9fJ] that, assuming statistical isotropy, 
these maps can be produced by convolving a normalized 
Gaussian white noise with the power spectrum of inter- 
est. The fact that we do not need to know anything else 



than the power spectrum of these effects should not be a 
surprise as the two-point correlation function encodes all 
the statistical information of Gaussian effects. Applying 
this technique with the spectra shown in the left panel 
of Fig. [SI we can therefore generate maps including the 
effects of all primary and secondary anisotropics. Now 
that we have these maps, we can use the second prop- 
erty mentioned above to include the effect of a network 
of cosmic strings. As the latter is uncorrelated with the 
primary and other secondary anisotropics, we indeed only 
have to add those maps to a map of the string-induced 
CMB anisotropics. The only unknown parameter con- 
trolling the overall amplitude of the string-induced tem- 
perature fluctuations is the string tension, for which we 
use GU = 7 x 10~ 7 (see Sec. HITB]) . We are then left 
with convolving the resulting map with the primary beam 
calculated in Sec. lIVBl While the tSZ will produce non- 
Gaussian spots, we can ignore them in this initial study. 
In practice, these sources can indeed easily be masked out 
of CMB maps when they correspond to bright (typically, 
5-c) sources. Moreover, in the case of multi-frequency 
experiments, it is possible to separate thermal SZ effects 
from thermal CMB distortions. 

We apply all these techniques to produce maps includ- 
ing the primary and all secondary anisotropics discussed 
above, as well as the effect of a network of cosmic strings. 
Figure [7] is given as a reference: it does not include any 
stringy effect. Therefore, the patterns associated with 
the acoustic peaks are clearly visible, as well as the Gaus- 
sian nature of the fluctuations. The right panel shows the 
corresponding normalized gradient magnitude, as defined 
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FIG. 8: CMB temperature anisotropics, including the same primary and secondary effects as in Fig. [7] as well as a string 
contribution of GU = 7 x 1(T 7 , for a 7.2° field of view at 147 GHz (left), and corresponding normalized temperature gradient 
magnitude map (right). Although the presence of strings can only be guessed in the left panel, very distinctive string features 
appear in the right one. In addition to sharp temperature steps scattered throughout the map, a few bright dipoles can be seen 
in localized string and loop regions with kinks and cusps. Once again, these maps take into account the effect of the primary 
ACT beam as derived in Sec. IIVBI 




FIG. 9: Normalized temperature gradient magnitude maps with the same components as in Fig. [8] but at 215 GHz (left) and 
279 GHz (right). The thermal Sunyaev-Zel'dovich effect vanishes at 215 GHz [see Eq. (|17[l ]. The enhancement of the string 
signatures with increasing frequency is clearly visible when comparing this figure to the right panel of Fig. [8] Moreover, when 
the latter is compared to the right panel of this figure, filaments produced by the tSZ become visible at 279 GHz. 
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in Eqs. (jT5j) and ([l6|l. with 9 Tes = 0.42'. By construction, 
the resulting map is highly non-Gaussian and enhances 
the rapid angular variations of the temperature patches 
of the fluctuation map. However, these structures re- 
main relatively smooth and, unlike the string patterns of 
Fig. do not exhibit any saturating bright spot. As a 
result, we were here able to use a linear color scale. 

In Fig. [5J we now include the effect of a network of 
strings with GU — 7 x 10~ 7 . The temperature fluctua- 
tions remain dominated by the Gaussian ACDM signal 
and the secondary anisotropies, with hints of tempera- 
ture steps induced by strings. Note the bright dipole in 
the bottom left corner of these maps: it is no longer 
symmetric in temperature due to the presence of the 
Gaussian sources. Even in the presence of the dominant 
primordial and secondary fluctuations, the gradient map 
still reveals the strings' location and velocity. However, 
the slow-moving segments are now hidden in the gradi- 
ent of the Gaussian patterns. As a result, the "faint" 
strings appear to be discontinuous. The normalized gra- 
dient color scale, although linear, has been limited to 30, 
while some of the saturating gradient values go up to 46, 
as one may have expected from the PDF shown in the 
left panel of Fig. [3] Finally, Fig. [9] shows the gradient 
magnitude of the temperature fluctuations generated by 
the same sources at 215 GHz (left panel) and 279 GHz 
(right panel), whereas the maps shown in Figs. [7] and [5] 
assume a frequency of 147 GHz. The temperature map 
is not represented since it is basically identical to the one 
of Fig. [5J The gradient color scale is kept the same as in 
this figure to make the comparison with the latter easier. 
As can be seen on these plots, the string "brightness" in- 
creases with increasing frequency, therefore reflecting the 
beam frequency dependence shown in the right panel of 
Fig.[5J As already mentioned, the increase in angular res- 
olution accompanying the increase in frequency enhances 
the gradient magnitude, but, interestingly, only for the 
intrinsically discontinuous anisotropies, such as the ones 
generated by strings. At these frequencies, the highest 
normalized gradient magnitudes observed go up to 65 at 
215 GHz and 95 for the 279 GHz channel. As can be seen 
by comparing the right panel of Fig. [5] to the left panel 
of Fig. [9l increasing the angular resolution also induces a 
higher sensitivity to the secondary sources. SZ features 
clearly appear on the 279 GHz gradient map as small fil- 
aments, whereas they are absent on the 147 GHz map, 
because of the extra smoothing coming from a narrower 
beam. The maps shown in this section therefore suggest 
that upcoming arcminute- resolution CMB experiments, 
such as ACT, could significantly improve the current con- 
straints on the string tension GU. In a subsequent paper, 
we will explore developing optimal non-Gaussian estima- 
tors aimed at detecting the coherent linclike string fea- 
tures in temperature and gradient maps. However, they 
are different enough from the features produced by the 
tSZ effect and other known sources of non-Gaussianity 
that we can already use a very naive statistic to get a 
crude estimate of the detection threshold that the afore- 



mentioned experiments could reach. Requiring that at 
least of order ten linelike features be visible in our gradi- 
ent maps, we can clearly see non-Gaussianitics down to 
GU ~ 2 x 1CU 7 . We suspect that more detailed studies 
will yield a firmer detection limit. 



V. CONCLUSION 

Based on numerical simulations of Nambu-Goto cosmic 
string networks, we produced 84 high-resolution string- 
induced CMB temperature maps over a 7.2° field. We 
extracted the expected probability distribution function 
and the angular power spectrum of string-induced CMB 
temperature anisotropies up to £ ~ 10 4 , including es- 
timates of the systematic and statistical errors associ- 
ated with our numerical simulations. Both of them ex- 
hibit distinctive features. In particular, strings should 
induce deviations from Gaussianity due to rare high tem- 
perature fluctuations, and lead to a power spectrum 
slowly decaying with increasing multipoles according to 
Eq. (TTJ} . Unless GU is well below the inflationary scale, 
the string-induced fluctuations start dominating the pri- 
mary anisotropies for £ > £ m - ln , with £ m - ln < 10 4 . How- 
ever, for values of GU smaller than a few times 10~ 7 , the 
string-induced anisotropies are themselves dominated by 
the tSZ, OV and nonlinear kSZ signals. 

We also discussed the observability of strings in a typi- 
cal arcminute-resolution experiment, focusing mainly on 
maps of the temperature gradient magnitude. Including 
the SZ effects in addition to the currently favored ACDM 
temperature anisotropies in our analysis, we find strings 
to be "eye visible," in the sense that at least of order ten 
string features are observable on a 7.2° gradient map, for 
GU down to 2 x 1CU 7 . As already mentioned, a drawback 
in our study is that we do not include the effect of point 
sources. The main risk associated with point sources 
could actually be in their removal from CMB maps. In- 
deed, for values of GU lower than 2 x 10~ 7 , the gradient 
maps no longer show the typical linelike shape of strings 
but still contain rare high temperature fluctuations com- 
ing from the active string regions. Since the latter are lo- 
calized, they may look like standard point sources, which 
could lead to their exclusion from CMB maps during the 
point-source removal process. A potential disambigua- 
tion method might rely on the intrinsic dipolar signature 
of these active string regions in the CMB temperature 
maps. Finally, let us mention that being able to take the 
gradient of a real temperature map requires a very good 
understanding of its noise properties, which could limit 
the applicability of the gradient technique described in 
this paper. 

In addition to the statistical analyses we discussed at 
the end of Sec. lTVCl other extensions of the present work 
might be worth thinking of. In particular, one might 
want to generate polarization maps, and especially their 
B-mode component since they have been shown to be 
string tracers [l6|, [9l| . Some of the results presented here 
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may also be relevant to the production of gravitational 
waves by cosmic strings. The presence of cusps is indeed 
expected to trigger gravitational wave bursts. Since they 
also lay on our past light cone, their number is certainly 
related to the number of bright dipoles one may observe 
in our simulated maps. 
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